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Spread of Infectious Diseases with a Latent 
Period 


Kanako Mizuno, Kazue Kudo 


Abstract Infectious diseases spread through human networks. Susceptible-Infected- 
Removed (SIR) model is one of the epidemic models to describe infection dynamics 
on a complex network connecting individuals. In the metapopulation SIR model, 
each node represents a population (group) which has many individuals. In this pa¬ 
per, we propose a modified metapopulation SIR model in which a latent period is 
taken into account. We call it SIIR model. We divide the infection period into two 
stages: an infected stage, which is the same as the previous model, and a seriously 
ill stage, in which individuals are infected and cannot move to the other popula¬ 
tions. The two infectious stages in our modified metapopulation SIR model produce 
a discontinuous final size distribution. Individuals in the infected stage spread the 
disease like individuals in the seriously ill stage and never recover directly, which 
makes an effective recovery rate smaller than the given recovery rate. 


1 Introduction 

Infectious diseases spread through human networks. Susceptible-Infected-Removed 
(SIR) model is one of the epidemic models to describe infection dynamics on a 
complex network connecting individuals. The ratio of the transmission rate to the 
recovery rate is called the basic reproduction number Rq. It is the expected number 
of infections caused by a typical infectious individual in a completely susceptible 
population |[Tl|2l. In the standard SIR model, the outbreak occurs when Rq > 1. 
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The likely magnitude of the outbreak, which is called the expected final size of the 
epidemic, depends only on Rq Elia. 

The spread of infectious diseases also depends on human mobility. In metapop¬ 
ulation SIR models, movements between different populations (groups) are taken 
into account iiiaia . Each node of the metapopulation network represents a group 
of individuals. Individuals can move between two nodes connected by a link. Al¬ 
though the epidemic threshold is Rq in each group, the global invasion threshold 
in the metapopulation system depends on the mobility rate as well as its network 
structure mM- 

In this paper, we propose a modified metapopulation SIR model in which a latent 
period is taken into account. Infected individuals behave like susceptible ones when 
they do not feel sick. They move between linked populations and spread diseases 
across different populations. We consider that such infected individuals are in a 
latent period. We assume that infected individuals become too sick to move after 
the latent period. Such ill individuals infect only the susceptible ones in the same 
population. This model is different from the SEIR model 13 , which is a common 
epidemic model in which a latent period is incorporated as an “Exposed” state. 
However, it belongs to a family of generalized SIR models that include multiple 
infectious stages El. The two infectious stages in our modified metapopulation SIR 
model produce a discontinuous final size distribution with a jump at Rq = 1. 

The rest of the paper is organized as follows. The metapopulation SIR model and 
the modified SIR model are introduced in Sec.|2] We demonstrate the discontinuous 
final size distribution of the modified model in Sec. |3] The effective recovery rate, 
which is different from the given recovery rate, is estimated, and it is the key to find 
what causes the discontinuity. Discussions and conclusions are given in Sec.|4] 


2 Model 

Eirst we introduce a metapopulation SIR model, which is an SIR model that is ex¬ 
tended to metapopulation networks. In the metapopulation SIR model, each node 
represents a population (group) which has many individuals, and each individual is 
in one of three states; S (susceptible), I (infected) or R (recovered). Individuals of 
state S are infected by those of state I in the same population. The infection rate is 
given by alm/Nm, where A,,, = Sm + hn-'r Rm with Sm, hn, and being the number 
of susceptible, infected, and recovered individuals of population m, respectively. In 
other words, the rate that S becomes I depends on the transmission rate a and the 
proportion of I in the same population. The constant rate that I becomes R, i.e., re¬ 
covery rate, is defined as j3. We here assume that all individuals move between the 
populations connected with links in the network at a constant rate w. The travel rate 
w is the same for all the individuals. The time evolution of the numbers of S, I and 
R in each population is described by 


Spread of Infectious Diseases with a Latent Period 3 

d,Sn = -aSnIn/N„+wY^{Sm-S„), (la) 

m 

dfln = CiSijIn/Nn jilnIn), (lb) 

m 

d,R„ = pIn + wY,{Rm-Rn), (Ic) 

m 


where the summations are taken over all the populations connected with population 
n. 

Next, we divide the infection period into two stages; an infected stage, which 
is the same as the previous model, and a seriously ill stage, in which individuals 
are infected and cannot move to the other populations. We call this model SIIR 
model. In this model, each individual is in one state of S (susceptible), H (infected), 
I (seriously ill), and R (recovered). Individuals of S in population m are infected and 
become H at rate a(Hm+Im)/Nm, where Nm = Sm+Hm+ht+Rm- Individuals of 
H become / at a constant rate /r. Individuals of I recover and become R at a rate j3. 
In the SIIR model, individuals of H move between the populations connected with 
links at a rate w, however, individuals of I do not. The time evolution of the numbers 
of S, H, I and R in each population is described by 


dfSn — —CXSnjHf, + In) /Nn + 'Wy'iSm — Sfi), ( 2 a) 

m 

^tHn — (^SniHn tn)/^n (2b) 

m 

dfln = ( 2 c) 

d,R„ = l5l„ + wY^{Rm-Rn), (2d) 

m 


where the summations are taken over all the populations connected with population 
n. 


3 Final Size Distribution 

The spread of a disease is expressed by attack ratio, which is the final proportion 
of R when I disappears in the entire metapopulation. The attack ratio plotted as 
the function of the basic reproduction number a/j3 is called a final size distribution. 
The final size distributions of the SIR model and SIIR models are shown in Fig. min 
this simulation, the number of individuals in each state is taken as a real number and 
the time step is discrete. We use a scale-free network with 900 nodes, whose degree 
distribution is P{k) ^ with 7 = 2.5. The essential results do not depend on 
gamma. In the initial state, 100 susceptible individuals belong to each node except 
for one randomly selected node in which one infected individual is included. The 
global invasion does not occur when a < j3 in the SIIR model as well as the SIR 
model. The change in attack ratio is continuous at a = j3 in the high-w region in the 


attack ratio 
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Fig. 1 Final size distributions of (a) metapopulation SIR model and (b) metapopulation SIIR model 
as the function of the transmission rate a and the travel rate w. In both cases, the recovery rate is 
j8 =0.5. 


SIR model, however, it is discontinuous in all region in the SIIR model. The shift of 
threshold in the low-w regions of the SIR model is often observed in metapopulation 
networks Eig. 

In this paper, we focus on the discontinuous final size distribution of the SIIR 
model. The jump in the attack ratio arises from the difference between the given 
recovery rate and an effective recovery rate. In the SIIR model, individuals H spread 
the disease like individuals I and never become R directly. Therefore, the effective 
recovery rate j3' is expected to be smaller than the given recovery rate j3. 

We show how to evaluate j3' below. Disregarding traveling between populations, 
the SIIR model (|2]i is rewritten as 


d,S = -aS{H + I), 
d,H = aS{H + I)-iiH, 
d,I = 
d,R = lil, 


(3a) 

(3b) 

(3c) 

(3d) 


where S = Sn/Nn, H ~ Hn/N„, I = In/Nn and R — R„/Nn. Combining Eqs. (iTbl i and 
(B, we have 

d,{H + I) = aS{H + I)-l5'{H + I), 




H + I 


P- 


We here take dfl = 0, which leads to // = (j3//i)/. Then, the effective recovery rate 
is calculated as 

Figure|2]illustrates that the evaluation of the effective recovery rate is appropriate. 
The simulation is performed in the same network with the same initial condition as 
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Fig. 2 The final size dis¬ 
tribution as the function of 
the transmission rate a for 
the SIR model with the given 
recovery rate j8 = 0.25 is com¬ 
pared with that for the SIIR 
model with the effective re¬ 
covery rate (3' = 0.25, which 
is calculated from Eq. 0 with 
j8 = 0.5 and fl = 0.5. Both 
curves agree in the region 
where a > 0.5. The travel rate 
w = 0.1 for both curves. 



Fig.IU The travel rate is w = 0.1, which is in the high-w region. The attack ratio 
for the SIIR model is calculated for j3 = 0.5 and ji = 0.5. In this case, the effective 
recovery rate is )3' = 0.25. The final size distribution for the SIR model with the 
given recovery rate j3 = 0.25 agrees with that for the SIIR model in the region 
where a > 0.5. This result implies the following. The effective recovery rate in the 
SIIR model is given by j3', however, global invasion cannot occur when a < j3. The 
difference between j3 and j5' causes the discontinuous final size distribution of the 
SIIR model. 

Since we disregarded traveling between populations when we evaluate the ef¬ 
fective recovery rate, the assumption that I is immobile should be irrelevant to the 
discontinuity in the final size distribution of the SIIR model. We now modify the 
SIIR model (l2]i, replacing Eq. (|2^ by 

d,In = IlHn - j3/„ + W^(/„ - /„). (5) 

m 

Figure|3]shows the final size distribution of the modified SIIR model. The simulation 
is performed in the same conditions as Fig.|2] The profile of the SIIR curve in Fig.|2] 
looks the same as the curve in Fig. [3 Therefore, the cause of the discontinuous 
final size distribution is the division of the infection period into two stages, and the 
mobility of I has no effect on the discontinuity. 


4 Discussions and Conclusions 


The effective recovery rate p', which is given by Eq. (|4]i, can be evaluated by another 
way. The basic reproduction number for the generalized SIR model that includes n 
infectious stages is given by 


^o = L 

/=l 


A’ 


(6) 
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Fig. 3 The final size distri¬ 
bution of the modified SIIR 
model in which H moves 
between populations, a is the 
transmission rate. The travel 
rate w = 0.1. 



a 


where a, is the transmission rate of the /th infectious stage, and 1 /jS,- is the mean 
duration of the stage miU. In our SIIR model, ai = a 2 = a, pi = ji and P 2 = P, 
and thus, Ro = a/ii + a/p ~ a{ii+p)/{iip). Therefore, 


n, a ^P 

^0 M + ^ 


(7) 


which is the same as Eq. dUl. 

In conclusion, the discontinuous final size distribution in the SIIR model is 
caused by the division of the infection period into two stages and the fact that the 
global invasion cannot occur when a < p. The final size distribution depends on the 
effective recovery rate p', and its shape coincides with that of the SIR model with a 
recovery rate p ~ P' in the region where a > p. 
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